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Liquid Front Profiles Affected by Entanglement-induced Slippage 
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Hydro dynamic slippage plays a crucial role in the flow dynamics of thin polymer films, as recently 
shown by the analysis of the profiles of liquid fronts. For long-chained polymer films it was reported 
that a deviation from a symmetric profile is a result of viscoelastic effects. In this Letter, however, 
evidence is given that merely a slip boundary condition at the solid/liquid interface can lead to an 
asymmetric profile. Dewetting experiments of entangled polymer melts on diverse substrates allow 
a direct comparison of rim morphologies. Variation of molecular weight M w clearly reveals that 
slippage increases dramatically above a certain M w and governs the shape of the rim. The results 
are in accordance with the theoretical description by de Gennes. 

PACS numbers: 68.15.+e, 83.50.Lh, 83.80.Sg, 47.61.-k 
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The controlled manipulation of small liquid volumes is 
one of the main tasks in the field of micro- and nanoflu- 
idics. In this context, the role of hydrodynamic slippage 
has attracted much attention due to its enormous rele- 
vance for applications such as lab-on-chip devices. A va- 
riety of experimental methods exist to characterize slip- 
page. The large number of techniques and results is sum- 
marized and discussed in recent review articles [H 0, 9 • 

If a thin liquid film on top of a solid substrate can gain 
energy by minimizing its interfacial area, the film rup- 
tures and holes occur (see [U, § and references therein). 
These dry circular patches grow with time, and the re- 
moved liquid accumulates in a rim surrounding the holes 
along their perimeter. Dewetting of thin liquid films on 
solid substrates is a result of internal capillary forces, 
which can be inferred from the effective interface poten- 
tial [H, Viscous dissipation within the liquid, which 
mainly occurs at the three-phase contact line y|, and 
friction at the solid/liquid interface (slippage) counter- 
act these driving forces. Concerning the characteristic 
dynamics of hole growth (hole radius R versus time t), 
Brochard-Wyart and colleagues Q proposed a sequence 
of different stages starting with the birth of a hole, which 
is characterized by an exponential increase of the dewet- 
ting velocity (denoted by V — dR/dt). It is followed by a 
regime that is governed by the formation of the rim and 
dominated by viscous flow (R oc t), followed by a stage 
of a "mature" rim, where surface tension rounds the rim 
and slippage might be involved. In case of full slippage, a 
growth law R oc i 2 / 3 is expected. This view was supple- 
mented by an experimental study by Damman and col- 
leagues, who associated different dewetting regimes with 
corresponding rim shapes in case of viscoelastic fluids @ . 
Vilmin etal. presented a theoretical approach ta king vis- 
coelastic properties and slippage into account In 
our earlier studies, we experimentally examined rim pro- 
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files of dewetting long-chained PS films on hydrophobized 
Si wafers and attributed the transition from symmetric 
("oscillatory") to increasingly asymmetric ("monotoni- 
cally decaying" ) shapes with the molecular weight of the 
polymer to viscoelastic effects [ll| . These results go along 
with a phcnomcnological model that predicts a phase di- 
agram for rim morphologies fl2l |. 

These previous studies, however, lack one facet, which 
by now has been identified as one of the major issues 
influencing dewetting dynamics and profiles: hydrody- 
namic slippage. In our recent studies we focused our 
interest on slippage of polymer melts, where polymers 
were used exclusively below their entanglement length 
0, G3, Ell- In this Letter, the studies are extended to 
long-chain polymer melts. Our aim is to characterize 
slippage of polystyrene films above their entanglement 
length. This might involve new problems related to non- 
Newtonian behavior due to viscoelastic flow and stress 
relaxation. The latter has been reported to cause poly- 
mer thin film rupture [l6l . [ItJ and to impact the early 
stage of dewetting [rjj. As will be shown in the follow- 
ing for the mature stage, residual stress and viscoelastic 
properties have no impact on the shape of dewetting rims, 
if low shear rates are involved. It will be demonstrated 
that in that stage, asymmetric rim profiles are caused by 
slippage only. 

Films of atactic polystyrene (PSS Mainz, Ger- 
many) with molecular weights ranging from 5.61 to 
390 kg/mol (termed e.g. PS(390k)) and low polydisper- 
sity (M w /M n — 1.02 — 1.09) were spin-cast from toluene 
solutions on freshly cleaved mica sheets and transferred 
to smooth hydrophobized Si wafers with a native oxide 
layer (obtained from Siltronic, Burghausen, Germany) by 
floating them on a Millipore™ water surface. To avoid 
residual stresses in films above the entanglement length 
of PS, the films were pre-annealed on the mica substrate 
well above their glass transition temperature (up to 3 
hours at 140 °C) QJ]. The prepared film thicknesses 
ho varied between 100 and 140 nm. Hydrophobization 
was either achieved through standard silanization tech- 
niques [HI using self-assembled monolayers of octadecyl- 
trichlorosilane (OTS) and dodecyltrichlorosilane (DTS) 



2 



TABLE I: Substrate properties. 



layer 


d (nm) 


rms (nm) 


dadv (°) 


AO (°) 


7,„(mW/m) 


AF 1600 


21(1) 


0.30(3) 


128(2) 


10 


15.0 


OTS 


2.3(2) 


0.09(1) 


116(1) 


6 


23.9 


DTS 


1.5(2) 


0.13(2) 


114(1) 


5 


26.4 



or by preparing a thin amorphous teflon layer (AF 1600) 
using the spin-coating technique. The thickness d and 
the root-mean-square (rms) roughness of the hydropho- 
bic layers (see Tab. QJ were measured via ellipsometry 
(EP 3 , Nanofilm, Gottingen, Germany) and atomic force 
microscopy (AFM, Multimode, Veeco Instruments, Santa 
Barbara, CA, USA). Additionally, substrates are char- 
acterized in terms of their wetting properties by their 
contact angles (advancing 9 a dv, and contact angle hys- 
teresis A9) of Millipore™water and their surface energy 
j sv (via measuring the contact angles of apolar liquids). 

Nucleation and subsequent growth of holes starts after 
heating the samples above the glass transition tempera- 
ture of the PS melt. In the following, optical microscopy 
is utilized to dynamically monitor the size of growing 
holes. At a certain hole radius (about 12 /xm) in the ma- 
ture regime, the experiment is stopped by quenching the 
sample to room temperature. The profile of the rim is 
then imaged by AFM [19| . As shown in Fig. [T]a), com- 
paring rim profiles for the same molecular weight, film 
thickness and viscosity (for the respective dewetting tem- 
perature) clearly exhibits strong differences that can ex- 
clusively be ascribed to the different substrates. For e.g. 
PS(65k) at 130 °C an oscillatory profile with a "trough" 
on the "wet" side of the rim is observed on the AF 1600 
substrate. A close-up of the trough is given in the in- 
set. Yet, in the case of OTS and DTS, rims are found 
that decay monotonically into the unperturbed polymer 
film. In Fig. [T]b), experimental profiles are shown for 
the AF 1600 substrate, where the molecular weight of the 
polymer melt is varied. In case of PS(65k), we observe 
a clearly oscillatory profile, whereas for e.g. PS(186k) 
or higher molecular weights, asymmetric (monotonically 
decaying) rim morphologies are recorded. Hence, an in- 
crease in M w on the identical substrate provokes the same 
morphological differences in rim shapes as before (Fig. [1] 

a) ) the variation of the substrate. To probe whether or 
not the effect is merely a result of the melt viscosity, 
the dewetting temperature on AF 1600 (see Tab. |TT| was 
additionally varied. Yet, no systematic influence of the 
viscosity on the rim profiles is found for this substrate, 
as exemplarily shown for PS(186k) in the inset of Fig.[T] 

b) . The observed marginal differences of the rim profiles 
can be related to small variations of accumulated mate- 
rial (via slightly different hole radii and film thickness in 
different experiments). 

Summarizing the qualitative observations, we find that 
an increase in molecular weight strongly amplifies the 
asymmetry of the profiles in the same manner as a 
change of substrate is capable of causing. According to 



a) 


600 




r 4 


135 
130 


1/ 


trni inh 

/ 9 

***** 




400 




/'•"*• \ 


125 


V 


Mr 1 OUU 


[nm 








3000 


4500 








\****« 








heigl 


200 




\ \ 
\ \ 














DTS 














OTS 















AF1600 


PS(65k) 


T=130°C 









2000 


4000 


6000 



distance \ from contact line [nm] 



b) 

600 





/~\ 400 




130'C 

/V. 140*C 

v — 150 " c 




/ *" N \ 

It v \ 200 

/' X \ 

A o 




PS(186k) on AF1600 




f •••X> ° 

\ " 


2000 4000 6000 




PS(271k) >v. 

-- - PS(186k) trou 9 n 

PS(65k) PSon AF1600, T=140°C 





2000 




4000 6000 



distance \ from contact line [nm] 

FIG. 1: a) Rim profiles of PS(65k) on DTS (dotted), OTS 
(dashed) and AF 1600 (solid line) at 130 °C. The inset illus- 
trates the oscillatory behavior on AF 1600. b) Rim shapes 
on AF1600 for different molecular weights (140 °C, PS(65k), 
PS(186k), PS(271k)). Inset: Variation of viscosity via differ- 
ent temperatures for the same molecular weight (see Tab. [TTJ) 
shows insignificant influence on rim profiles. 



our previous studies on OTS and DTS concerning non- 
entangled PS melts, asymmetric rims are related to slip- 
page [ll, [3, [H. For DTS, slip lengths b of up to 5 
/jm were obtained, for OTS, values are about one or- 
der of magnitude smaller. The question now arises as to 
the manner in which slippage correlates to the molecular 
weight and the number of entanglements in the AF 1600 
system. To quantify slippage there, rim profiles will be 
analyzed. 

As shown recently, lubrication models for Newtonian 
flow and different slip conditions (no-slip, weak-slip, 
intermediate-slip and strong-slip [2(|) and a more gener- 
alized model based on the full Stokes equations [14J can 
be applied to thin dewetting films with the important 
result that values for the slip length b can be obtained. 

To assess whether these models might be used to an- 
alyze the experimental data, the impact of viscoelas- 
ticity on fluid flow needs to be estimated: The Weis- 
senberg number Wi = rj relates the time scale of stress 
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TABLE II: Comparison of viscosities (obtained from the WLF 
equation), calculated relaxation times r and measured rim 
widths w (on AF 1600) given for a selection (see Fig. [1} of 
dewetting temperatures T and molecular weights M w . 



M w 


T (°C) 


•q (Pa s) 


r (s) 


w (/^m) 


65k 


130 


1.3 x 10" 


6.5 


2.73 


65k 


140 


1.6 x 10 5 


0.8 


2.81 


186k 


130 


4.4 x 10 7 


222 


5.02 


186k 


140 


5.6 x 10 6 


28.1 


4.86 


186k 


150 


1.1 x 10 6 


5.4 


5.29 


271k 


140 


2.0 x 10 7 


100.5 


7.47 



relaxation r with the time scale of applied shear 7. 
Relaxation times are determined from the correspond- 
ing viscosities r\ and the shear modulus G rj 0.2 MPa 
(for PS, HU) via r = 77/G (see Tab. HTJ) . Evaluation 
of the shear rates within the moving rim is based on 
7 = max(d x u x , d z u x ) « max(s/Ax,s/Az), where s is 
the dewetting speed (identified as the dominant velocity 
u x 3> u z within the rim) which can be inferred from opti- 
cal hole radius measurements. The lateral displacement 
Ax is estimated as the width w of the rim. The value w is 
defined as the distance between the three-phase contact 
line and the position where the rim height has dropped 
to 110% of the film thickness ho, i.e. h(w) = l.l/io- The 
vertical displacement Az is given by Az = ho + b, where b 
denotes the slip length as determined by the quantitative 
rim analysis described in the last part of this Letter. The 
described estimation leads to values for Wi between 0.01 
and 0.1 (on AF 1600) and 0.5 (on OTS/DTS). An influ- 
ence of viscoelastic effects on flow dynamics is expected 
for Weissenberg numbers of 1 and larger. In that case, a 
model including viscoelasticity has to be used [22]. Low 
shear rates as described above and the pre-annealing step 
after spin-coating allow us to safely exclude viscoelastic 
effects as the source of different rim morphologies. This 
consideration justifies the assumption of Newtonian flow, 
on which the following data analysis is based. 

The theoretical description of the model for Newtonian 
liquids presented in [14j is derived from the full Stokes 
equations in two dimensions for a viscous, incompressible 
liquid including the Navier slip-boundary condition b — 
u/d z u\ z =o, where u denotes the velocity in x-direction 
and z the vertical scale. This leads to the equations of 
motion for a flat liquid film. Then, a linear stability 
analysis is applied by introducing a small perturbation 
Sh = h(x, t) — hg and velocity u{x, t) to the undisturbed 
state h — ho and u = 0. In the co-moving frame £ = 
x — s(t), where s(t) is the position of the rim, a quasi- 
stationary profile evolves. Solving the linear equation 
using the normal modes ansatz 5h — Shoexp(k£ i ) and 
u = uq exp(fc£) gives a characteristic equation for k which 
is expanded up to third order according to the Taylor 
formalism: 
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FIG. 2: Slip lengths of PS films on AF 1600 obtained from 
rim profile analysis. Each data point represents the average 
value of different measurements (2-3 measurements for each 
M w value above 35.6k and 5 or more below). The solid line 
represents a fit to the experimental data from PS(35.6k) to 
PS (390k), the dashed line visualizes b — aN 3 /N^ using values 
for a and N e given in the literature The dotted line 

represents the average of slip lengths for melts from PS(5.61k) 
to PS(18k). 



(1 + ||)(^fc) 3 + 4<7a(l + ^)(hok) 2 - Ca^ = 0, (1) 

where Ca denotes the capillary number given by Ca = 
V^Msi- Hence, rim profiles depend on the capillary num- 
ber Ca and the ratio of slip length b to film thickness 
ho. Eq. ([1]) implies that a morphological transition from 
oscillatory (complex conjugate solutions for k) to mono- 
tonically (real k values) decaying rim shapes occurs if 
Ca 2 > (3 3 (b/h + l/3) 2 )/(4 4 (b/h + 1/2) 3 ) on slippery 
substrates (large b/ho) and for large capillary numbers 
Ca (l4lll5j|. Furthermore, Eq. ([I]) predicts progressively 
asymmetric rims for increasing values of b/ho. 

In the following, PS melts on the AF 1600 substrate 
are considered: Rim profiles are evaluated by fitting the 
wet side of the rim with a damped oscillation (oscillatory 
profiles, up to PS(81k)) or a simple exponential decay 
(monotonic rims, for PS(101k) and above). 

In the case of an exponential decay, a single value for 
k is gained from the experimental data. Knowing the 
film thickness ho and the capillary number Ca (from in- 
dependent rim velocity data), this k value allows us to 
calculate the respective slip length via Eq. ((T|). For oscil- 
latory profiles, k is a pair of complex conjugate numbers. 
From these two independent solutions, both the capil- 
lary number and the slip length can be extracted via Eq. 
([T]) at the same time. Viscosity data calculated from 
these Ca values showed excellent agreement with theo- 
retical values (from the WLF equation) and demonstrate 
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the consistency of the applied model. For further details 
concerning the fitting procedure and data evaluation, one 
should refer to previous publications [HI, 0, EH] ■ 

The slip lengths 6 of a set of experiments, i.e. differ- 
ent dewetting temperatures and viscosities, are shown in 
Fig. Has a function of molecular weight M w of the PS 
melt. Two distinct regimes can clearly be seen: Con- 
stant slip lengths in the order of 10 to 100 nm at small 
M w and slip lengths that increase with the power of 
2.9(2) in case of larger M w . The transition between both 
regimes occurs at 37kg/mol, which is consistent with the 
critical value M c of entanglement effects (M c =35 kg/mol 
for PS according to [2l[). For the same viscosity (e.g. 
2 x 10 6 Pas) in case of e.g. PS(51.5k) and PS(271k), 
an enormous difference in slippage of two orders of mag- 
nitude solely due to the difference in chain length is ob- 
tained. We therefore state that slippage is directly re- 
lated to entanglements rather than viscosity. Besides the 
fact that an influence of the melt viscosity (via temper- 
ature variation) for a certain molecular weight could not 
be detected, even experiments omitting the pre-annealing 
step after preparation on mica do not show significant 
impact. Additionally, we want to emphasize that careful 
analysis of the dewetting velocity data from hole growth 
dynamics (of the identical samples), according to the 
model presented in (23j, confirmed the b(M w ) behavior 
obtained from rim analysis [24j . 

For long polymer chains, described by the reptation 
model, de Gennes predicts a scaling of b oc iV 3 , where 
N denotes the number of monomers, which corroborates 
the experimental exponent of 2.9(2) for experiments with 
M w > 35 kg/mol as depicted in Fig. [2j According to the 
model, b = aN 3 /N^, where N e is the bulk entanglement 
length (number of monomers in an entanglement strand) 
and a a polymer specific molecular size, and is plotted 
in Fig. [2] as the dashed line for the literature values of 



= 163 [21] and a = ZA |H for PS. From the lin- 
ear fit to our experimental data (solid line in Fig. 
we obtain a/N^ =1.12(7) x 10 _5 A, which is one order of 
magnitude smaller that expected by using the literature 
values given above. Keeping a from the literature, a sig- 
nificantly larger entanglement length of N e = 517 results 
from the linear fit. This is an indication for a significantly 
lowered effective entanglement density in the vicinity of 
the substrate compared to the bulk and is in line with 
recent studies of thin polymer films [13, HH . According 
to Brown and Russell, N e of a polymer melt near an in- 
terface is expected to be about 4 times that in the bulk 
[29l ] . which corroborates our experimental results. 

To conclude, entanglements can strongly amplify slip- 
page as deduced from dewetting experiments. Assigning 
asymmetric profiles exclusively to non-Newtonian effects 
such as viscoelasticity therefore turns out to be invalid. 
On the contrary, the hydrodynamic boundary condition, 
i.e. slippage, can dramatically provoke morphological 
changes in rim shape such as transitions from oscillatory 
to monotonic profiles. The onset of slippage correlates to 
the critical chain length for entanglements. The depen- 
dence of the slip length on molecular weight corroborates 
the description by de Gennes. To derive vice versa slip 
lengths solely from de Gennes' theory entails that bulk 
values for a and N e (from literature) have to be used that 
might not appropriately describe polymer melts next to 
an interface. The method described here, however, rep- 
resents a powerful tool to gain experimentally access to 
liquid flow properties at the solid/liquid interface and to 
quantify slip lengths. 
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